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Abstract 

We characterize the asymptotic performance of a class of positive operator valued measurements 
(POVMs) where the only task is to make measurements on independent and identically distributed 
quantum states on finite-dimensional systems. The POVMs we utilize here can be efficiently described 
in terms of a reasonably small set of parameters. Their analysis furthers the development of a quantum 
method of types. They deliver provably optimal performance in asymmetric hypothesis testing and 
in the transmission of classical messages over quantum channels. 

We now relate them to the recently developed a — z divergences D a ,z by giving an operational 
interpretation for the limiting case liiriQ.-ji D a ^- a in terms of probabilities for certain measurement 
outcomes. This explains one of the more surprising findings of [Tj in terms of the theory of group 
representations. In addition, we provide a Cauchy-Binet type formula for unitary matrices which 
connects the underlying representation theoretic objects to partial sums of the entries of unitary 
matrices. 

At last, we concentrate on the special case of qubits. We are able to give a complete description of the 
asymptotic detection probabilities for all POVM elements described here. We take the opportunity to 
define a family of functions on pairs of semi-definite matrices which obeys the quantum generalizations 
of Renyi’s axioms except from the generalized mean value axiom. This family is described by limiting 
values of a — z divergences for the extremal values of the parameter. 


I Introduction 

The importance of representation theory for quantum information is understood best by taking a quick 
look at the structure of communication systems: Throughout, these systems employ certain structures 
which are inserted into them by construction at sender’s side and can then be detected by the receiver 
even if the signal gets corrupted by noise. The most general approach for making a signal received 
despite noise is to use some form of repetition. Signals arising from repetition are obviously invariant 
under permutations. In a probabilistic sense, this property may continue to hold even under the influence 
of noise. Subsequently, the early and simple intuition of making communication resilient against noise 
by repetition has been developed to what is modern communication theory. The simple invariance under 
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permutations that repetition delivers got lost in the hunt for higher throughput but the importance of 
group actions in code design prevailed. 

Of course this argument is independent of the mathematical model we use for our signals and so an 
analogous reasoning applies to the study of signals which are described by quantum theory. 

We concentrate on the study of finite dimensional quantum systems here. Measurements on such systems 
are completely described by positive operator valued measurements (POVMs). Such measurements have 
a finite number of different possible measurement outcomes. The probability that a POVM yields a 
specific outcome depends on the state that the system is in. In order to make communication tasks 
viable it is of great importance to deliver both sets {p m } mE m of signals which embody the messages 
to G M into the quantum states p m and detection schemes {D m } me m in the sense of POVMs such that 
the probability p m of getting measurement outcome to when the signal state is to satisfies p m ss 1. 
Besides these requirements it is necessary to deliver efficient descriptions of both the signals and the 
measurements. 

We take this as an easy to grasp motivation to study the asymptotic performance of certain types of 
POVMs which are built up from certain well-described representations. This provides a way to obtain a 
clear specification of a class of quantum measurements in terms of a reasonably small (as compared to the 
frequency typical subsets of classical information theory. Our description needs less than twice the number 
of parameters needed for the description of a frequency typical subset with the underlying alphabet being 
[d]) set of parameters, an approach which opens up the possibility to deliver standardized detection 
procedures for quantum communication. Our main contribution in this direction is the characterization 
of the asymptotic detection probabilities of the POVMs studied here in terms of convex optimization 
problems which require optimization only on a finite number of copies of the underlying systems. 

The POVM elements that we analyse here are capable of delivering provably optimal results both in 
asymmetric hypothesis testing |20| and in message transmission over noisy channels [5]. They also deliver 
an intuitively appealing step forward in the development of a “quantum method of types”. In this work we 
are further able to show initial connections to the geometry of positive matrices and to prove a surprising 
connection to the recently developed a — z divergences [I] . 

In the work (20] the author defined this specific class of POVMs as follows: For a given orthonormal 
basis {ej}f =1 of C d and a “frequency” or “type” (a nonnegative function / : {1,..., d} —> N satisfying 
X], /(*) = n), consider the irreducible representations of the symmetric group S n on the frequency typical 
subspaces 

Vf := span({eq (8)... ® d n : |{/c : ik = i}\ = f(i) for all i}. 

Since these subspaces are invariant under permutations by definition, they naturally split up into different 
isotypical subrepresentations Vf t x (where A denote Young frames and some Vf t x nray not contribute to 
above decomposition, meaning that Vf t x = {0} for some pairs (/, A), while some representations may 
occur several times) 

Vf = ®V f , x. ( 1 ) 

A 

The reason that the projections Pf,x onto these subspaces deliver optimal results in asymmetric hypothesis 
testing stems from the following observation: Given a state a with eigenvalues ti > t^ > ... > td, > 0, we 
may pick one of its eigenbases for the definition of the Vf. It is then straightforward to show that the 
estimate 

ti{P f , x a® n } » dim(F A ) • 2~ n ^ TifWogU 

is valid, where dim(i 7 A ) is the dimension of the irreducible representation F x of S n corresponding to A 
and H(iA) = — (T),: -Ajlog(-Ai) is the entropy of the normalized Young frame. It turns out that, for 
A « n • spec(p) and f(i) ~ n ■ (e*, pet) V i € [d] for some arbitrary second state p we get 

tT{P fi xa® n }^2- nD< ' p ^'>. (4) 


( 2 ) 

(3) 
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In our earlier work [20j we have not been able to characterize the exact asymptotic behaviour of the maps 
(/, A) !->■ trjP/^cr®"} when the basis {ei,..., ej} is no longer such that the representation of cr in that 
basis is diagonal. We now start investigating this topic: For a fixed basis, two probability distributions 
p and s on {1,..., d} for which s(l) > ... > s(d ) we are interested in the map 

(cr,p, s)h>- lim i logtr{P /( „) A( „)cr® n } (5) 

n—> oo 77, 

when —> p and —> s. The above limiting procedure is unfortunately not always easy to 

characterize. We do therefore not give a characterization in all detail and rather concentrate on the cases 
where /( n ' ) = Afor all n € N or where a is diagonal in the chosen basis, whenever d > 2. Only in 
the case d = 2 we are able to exploit the specific nature of qubit systems to deliver a more detailed 
description. 

Since every state a can be transformed into a state UaU ^ which is diagonal in {ei,..., ed} we get a 
representation of a as 

a ®u = ^ ( 6 ) 

f \ ' f' 


even if a is not diagonal in {ei,..., e^}. Thus we felt motivated to also study the asymptotic behaviour 
of Hilbert-Schmidt scalar products of Pf t \ with A® n Pp x A'® n for arbitrary operators A. Especially in 
the case where A = U is a unitary this led us to prove an interesting algebraic formula (see Lemma [3]) 
which arises from the study of the minimum of the function 


( U,p,q,s ) 


lim — 

n —¥oo 77, 


log tr{P /( „) iA(B) U® n P g(n) iA(B) P t0n } 


(7) 


where -f^ n > —> p, A g 1") —> q and Aaw —> s. A detailed description of our approach is postponed to 
Section ED where we also give precise definitions of our main objects. 

The method we utilize here opens up the possibility to split the analysis of the detection procedure 
into two parts: We note that [P/,P\] = 0 for all / and A. Thus, one can always realize e.g. the P/ 
measurement first. This task is comparable to detection procedures in classical systems, if one takes local 
measurements (e.g. measurements where each POVM element is of the form ®™ =1 M; for Mi,..., M n > 0 
and Mi,..., M n £ Md for free. It is not as clear how to implement the P\ measurements. A method for 
doing so has been outlined in m, with its success being conditioned on the physical realizability of what 
is called the Schur-transform there (and in [2]). Apart from [TTj . the work [B] gives a lot of structural 
insights into the relations between representation theory and quantum information theory. 

Recent work has also put into focus the definition of quantum relative entropies, and large families of 
such quantities have been defined and key properties like unitary invariance, convexity or monotonicity 
have been proven to hold. 

We do not make any attempt to give a complete overview on the topic, we rather point the reader to the 
papers mmm- These contain a good amount of the necessary history as well. Our focus here will be 
on the notion of a — z relative entropies D az which are defined for real parameters a/1 and z ^ 0 in 
PQ. It was proven by the authors of that work that their definition includes all the previous ones in the 
sense that certain choices of the parameters a and 2 yield the other relative entropies. 

We will pay special attention to the limits linia^i D 0} i_ a and limo^i D a z for z ^ 0 which yield a rather 
cumbersome formula in the first case and the quantum relative entropy in the second. We are able here 
to give a direct operational interpretation of the first quantity in terms of an asymptotic probability (for 
details, see Section m of obtaining certain measurement outcomes when the POVM which is utilized 
is of the form {Pf,x}f,x for some choice of basis. The second quantity is already connected via |20j and 
[2Tj . This motivates our definition of a huge set of non-negative functions R. which fulfill some of the 
Renyi axioms (but not the generalized mean value axiom) and can be parameterized in such a way that 
they naturally include both lim Q _»i £) aj i_ a and limo^i D a z for z / 0. We point out a possible way to 
derive further connections to the geometry of positive matrices in Section [Vi] 

A further study of this interplay between representation theory, information theory and the geometry of 
positive matrices is postponed to future work. 
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II Notation 


Throughout, d £ N denotes the dimension of the vector space C d that we will be working on. We use the 
convention [d] := (1,..., d}. For two natural numbers d and k satisfying k < d and every j £ [d] we set 
[d\j := {(«i,... ,ifc) : 3! I £ [d] : ii = j}. The set of probability distributions on [d] is *}3([d]), and the set 
of ordered elements of fp([d]) is *}}l([d]) := {p £ f)3([d]) : p( 1) > ... > p(d)}. For an arbitrary p £ fp([d]), 
pi g Cpf-Qd]) is defined to have the same values as p, but in descending order. A function g : [d] —>■ R 
satisfies g h f for another function / : [d] —► R if J2i=i s(*) — Yli= i /(*) for all k = 1,..., d. 

The set of positive matrices acting on C d is Vd , the set of matrices is M,j. Non-square matrices are 
elements of Mdxd '> where d is the number of rows and d' the number of columns. The adjoint of b £ Md 
written b'. 

S(C d ) is the set of states, i.e. positive semi-definite matrices with trace (the trace function on Md 
is denoted as tr) 1 acting on the Hilbert space C d . Pure states are given by projections onto one- 
dimensional subspaces. A vector x £ C d of length one spanning such a subspace will therefore be referred 
to as a state vector, the corresponding state will be written |x)(x|. For a finite set X, |X| denotes its 
cardinality. If X' C X, then X\X' := {x £ X : x ^ X'}. The extremal points of the convex set iP(X) 
are the distributions S x defined by S x (x') = 1 if and only if x = x’. We will also need “channels”, 
meaning probability preserving linear maps from *jI(X) to ip(Y). These are represented by matrices 
W = (■u;(y|a:)) x6 x,yeY which satisfy Y1 V £y w (y\ x ) = 1 f° r & H £ £ X. Their action is uniquely defined by 
setting W(S X ) := w(y\x)S y for every a; £ X. The set of channels from X to Y is denoted C(X, Y).j 

For any n £ N, we define X” := {(xi,..., x n ) : x* £ X Vi £ {1,..., n}}, we also write x n for the elements 
of X n . Given such element, N(-\x n ) denotes its type, and is defined through Vx £ X : N(x\x n ) := |{i : 
Xi = x}\. The set of all types arising from words of length n is written T„ or, if the alphabet is not clear 
from the context, T„(X). 

The von Neumann entropy of a state p £ S(H) is given by 

S(p) ■= -tr(plogp), (8) 

where log(-) denotes the base two logarithm which is used throughout the paper. The entropy of r £ fP(X) 
is defined by the formula 


H ( r ) ■= - X] r{x)log(r(x)). 

aiGX 


(9) 


Given two states p, cr £ S(C d ), the relative entropy of them is defined as 


D(p\\cr) 


tr{p(log(p) - log(cr))}, if supp(p) C supp(cr), 
oo, else 


( 10 ) 


For p, q £ fP([c?]) we may set p := X)f=i P(®)l e *)( e *l an d a := Z)f=i <7(*)l e *)( e *li then (with a slight abuse of 
notation ) D(p\\q) := D(p\\a) defines the classical Kullback-Leibler distance D(p\\q) between probability 
distributions p,q £ ^([d]) as well. 

We now fix our notation for representation theoretic objects and state some basic facts. 

The symbols A, p will be used to denote Young frames. The set of Young frames with at most d £ N rows 
and n £ N boxes is denoted Y d, n - 

For any given n, the representation of S n we will consider is the standard representation on (C d )®" that 
acts by permuting tensor factors. 

The most important technical definition for this work is that of frequency-typical subspaces Vf of (C d )® n . 
These arise from choosing a fixed orthonormal basis {e^}^ of C d , choosing a frequency / (a function 
/ : [d] -t N satisfying Yh=i /(*) = n )> setting T f := {(i 1? ... ,i n ) : |{4 : i k = j}| = /(j) Vj £ [d]}, and 
defining 


Vf := spande^ <8»... <g) e in : (*i, ...,*„) £ Tf}). 


( 11 ) 
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They have been widely used in quantum information theory, but share one very nice property that 
has not been explicitly exploited in quantum information theory until BUI : They are invariant under 
permutations, if the (linear) action B of S n on (C d )® n is defined in the natural way via 

B(r )v\ ® ® v n ■— iv-i(i) ® ® t» T -i(„) (12) 

for all t £ S n and Vi,... ,v n £ C d . From the invariance of each Vf under the action B of S n it immediately 
follows that 


V/ = ®V f , a , (13) 

A 

where each V/ t \ is just a direct sum of irreducible representations corresponding to A that is contained 
entirely within Vf. The multiplicity of F\ within Vf is given by dim(Vy >- \)/dim(F\). It is a number 
which scales at most polynomially in n, if d is kept fixed. The quantity F\ denotes the unique complex 
vector space carrying the irreducible representation of S n corresponding to a Young frame A. Each such 
A consists of n £ N boxes and has row lengths Ai,..., A d for some d £ N. Thus, A defined by A(*) := A A,; 
for every i £ [d] defines an element of *}jl([d]). 


During our analysis it turns out that, for every k < d, the vectors 


Vk ■= 


- 7 = ^2 sgn(r)B(r)ei <g) ... <g> e k 
V k res k 


(14) 


are important. Also, we are going to employ the following estimate taken from [Tj Lennna 2.3], which is 
valid for all frequencies / : [d] —► N that satisfy J2i=i /(*) = n: 


1 

(n + l) d 


< \T f \ < 2 nHU \ 


(15) 


where / := i/. We will also need [7| ( Lemma 2.7] which employs the variational distance that we define 
as ||p — <?|| := \p( x ) — d(®)| fo r all p, q £ ^3(X) and delivers: 

Lemma 1. If, for X a finite alphabet and p,q £ ')l(X) we have ||p — <?|| < 0 < 1/2, then 


\H(p)-H(q) |<-01ogj||. (16) 

Another very important estimate is the following one (a derivation can e.g. be found in ]I9]): 

2 n(H(A)-^log(2„)) < dimFA < 2 nH( A) ( A £ y djn ). (17) 


During our investigation we shall need the following sets of distributions: For every q £ f)3([d]) and 
1 < k < d, set 




p£¥({d] k ) 


N(i\(i i,.. .,i k )) > 1 => ■ ■ ■,ik )) = 0 

P([d]i) = k ■ q(i) V i £ [d] 


(18) 


Such distributions can be constructed by taking a unitary U £ Md and defining p £ < )3([d] fe ) via 
p((*i,..., ik)) ■= |(<g> ... ei k , U® k v k )\ 2 . This ensures the validity of ... , i k )) >1 => 

p((*i,..., ik)) = 0. Lemma [3] then delivers the values of p([d}^). This connection demonstrates that 
^(q, k) ^ 0 is possible, thus making our definition nontrivial. 

We now switch the topic one last time in this section and concentrate on additional entropic quantities 
which are necessary in the remainder: 
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Definition 1 (Reverse Sandwiched Relative Entropy). For p,cr £ S(C d ) with suppp C supper and 
a £ R\{0}, set 

D a (p\\cr) := j- log tr ^ ap 2 ^ a - a '> ^ |. (19) 

The revere sandwiched relative entropy is derived from the sandwiched relative entropy, which was 
defined in [23] and [13] : 

Definition 2 (Sandwiched Relative Entropy). For p, a £ S(C d ) with suppp C supper and a £ R\{0}, 
set 

D a {p\W) ■= ^-j-logtr | (a^pa^^j |. (20) 

As was made explicit in [I] Equation 11], the two quantities are related through the following equation: 

(a - l)D a (p\\<j) = (—a)Di- a {(j\\p). (21) 

The sandwiched relative entropy D a has been proven to have a huge number of highly desirable properties 
in, among others, the work [3, and [18] . From one of its origins, it is intimately connected to quantum 
channel coding [23]. Applications are also found in hypothesis testing EE], ini, mi, m\- 
A more general definition was made by Audenaert and Datta |li. It adds the parameter z > 0 and reads 

Definition 3 (a — z relative entropy). For p,a £ S(C d ) with suppp C supper, a £ K\{0} and z > 0, set 

d *Ap\W) := ^rylogtr {(p a/z v (1 - a)/z ) Z }. (22) 

The work jT. not only defined this quantity but also provided a lot of details and especially gave an 
explicit formula for the limit lim a _»i _D aj i_ a (/?||cr) that we shall use in the proof of Theorem^ 

III Main Results and Definitions 

Throughout, we make our definitions with respect to one fixed but arbitrary orthonormal basis Bd = 
{ei,..., edj - the standard basis within C d . Every matrix and also every quantum state are represented 
with respect to that basis. 

Definition 4. To any pair ( p , a) £ S(C d ) xS(C d ) we assign a unitary transformation U p such that U p pUJ 
is diagonal and has its diagonal entries sorted in descending order. We can then define the function 

4> : S{C d ) x S(C d ) -)• R+, (23) 

(p, cr) h-> - lim - logtr {U}® n P fM>xM U® n CT® n }, (24) 

n—>• OO n 

where the sequence (/^ n \ A^ n ^) ra£ N satisfies linin^oo A-A*-") = ((ei,pei),..., (e^, ped) and /*■") = A^ for 
all n £ N. 

Although this definition is ambiguous whenever p has degenerate eigenvalues we show later that it is 
still well-defined. 

Definition 5. In the same way as in Definition we take any unitary matrix U a such that U^aUa is 
diagonal in Bd- We define the function A : <S(C d ) x S(C d ) —> R+ by taking any sequence of frequencies 
satisfying lim,,-^ A/( n ) = ((ei, U a pU\ei ),..., (e^, U a pU\ed)) and a sequence of Young frames satisfying 
liinn^oo tJ-A^ = spec(p). We then set 

A : S(C d ) x <S(C d ) —>• M+, (25) 

(p,a) ^ - lim - logtr {U® n P fM xM U\® n a® n }. (26) 

n —>oo 77, J 

Again, we show later that this definition does not depend on a particular one among the many possible 
choices of U a • 
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The third quantity we define is 


Definition 6. Let d £ N and A £ Md. Define a function as 

A a : *p([d]) x ^{[d]) x S(C d ) -)■ R+, 


(27) 


- lin Woo i log tv{A® n P fM A( „) if s >- 

oo, else (28) 

where (/ (n) )„ 6 n and (A (n )) neN are sequences satisfying lim^oo = s £ ^([d]), lim n _ ! . 00 i/ (n) = 

p £ tp([d]) and the sequences are constructed such that A ^ A for all n £ N. 

Of course A u = $ whenever U is a unitary matrix and U pU' is diagonal in Bd and has decreasing 
diagonal entries and s = spec p. Also, A jj = A whenever U is unitary and U'crU is diagonal in B and 
p = pinch p, s = spec p. In the remaining cases it is not clear from their definition that <f>, A or A are 
well defined. Note that implies the existence of sequences (A^ n )) n6 N and /^) n e n with respective 

limits s and p and such that the Kostka numbers K X ( n ) jWi- of these sequences are non-negative [TO. 
Exercise 2]). If K X ( n )j(n > > 0 however then the construction provided in [22] Chapter 5.5] proves that 
Vf(n) X in) ^ {0}. We will use this connection more explicitly in the proofs of Theorems Q] and [4] 

The question whether A a is well-defined in the sense of being independent from the specific sequences 
(/W) n6 n and (A<">) n6N will be settled here only for d = 2 in Theorem 0] We are thus left with a 
conjecture: 

Conjecture 1. The functions Aa from Definition 0 are well-defined for every d > 1 and every A £ Md- 

For the other two quantities it will become immediate that they are well-defined once we calculate 
the limits in the Definitions [4] and [5] This task leads us to the following theorem: 

Theorem 1. For every two states p,a £ S(C d ) we have 

1. $(p|k) = lim D a {p\\a), 


(P, 8, cr) >-)■ 


2. A(p\\a) = D{p\\a). 

Remark 1. The second of the above statements has been proven in f20 )/ and will not be proven here 
again. 

This result raises some interest into a more in-depth study of the projections P/ jA . As we already 
observed before [20] the subspaces V/ )A are generically not irreducible. In these cases computations are 
less straightforward as in the cases where we have irreducible representations. We therefore concentrate 
here on a study of cases where at least one of the subspaces involved into the calculation is irreducible. 
The decomposition 


= V tr{P f (UaU^ n } 
^ IT, I 


x jj® n 


(29) 


which is valid whenever cr is diagonal in Bd motivates the study of objects of the form 
tr{Pf )X U® n Pf^ X 'U^ n } = tr{Pf tX U® n Pf', x U*® n }6{X, A') where U £ M d is unitary. Of course this is 
done again in the asymptotic setting, and with a slight increase in generality: 

Definition 7. For every d> 2, asymptotic shapes s of Young frame, frequencies q £ ?P([d]), and matrix 
A £ Md we define 

Q(q,s,A) := - lim -logtr{P A( „> xM A® n P f(m xM A^ n }, (30) 

n —Too n 

where (X^) n eN is a sequence of Young frames satisfying hin^^oo j^X^ = s and (/^) n gN a sequence of 
frequencies satisfying lim^-yoo X f (") = g. 
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The results we obtain from the study of 0 are presented in the next theorem: 

Theorem 2. For every d > 2 and asymptotic shape s of Young frames as well as asymptotic frequency 
q and A G Md the function 0 assumes the value 


d 

Q(q, s, A) = H(s) — min V'' s(k) min \D(p\\p k ) (31) 

1 ; w W(s)= q f^ y j P GV(W(S k ),k) k K ’ 

where p k (h,.. .,i k ) := \{e ix <g).. -®e ik , J2 T eS k sgn(r)B(r) (ei <g) .. . e fe )}| 2 and s G ^([d]) is defined 

by s(k ) := (s(k) — s(k + 1)) • k for all k G [d\ using the convention s(d+ 1) := 0. For every unitary matrix 
U and fixed asymptotic shape s, the function q K► Q(q, s, U) assumes its minimum at a distribution q 

which satisfies q(i) := Ylt-i s(k) i ^ U ' k ■ 

Al I 12 

Remark 2. Note that the map k H > i~~ actually defines an element of C([d], [d]) since for every 

k we have E?=i Ji T L = 1 - 

Especially the location of the minimum which we describe above made us conjecture an interesting 
formula via the following route: An application of Pinsker’s inequality to above formula for 0 lets us 
transform the search for the minimum into a question about distance in norm rather than relative entropy. 
This decomposition delivers a lower bound which can be shown to equal zero if and only if q has the 
desired form: 

Lemma 2 (Estimate for norms). Let s G ^([d]) and define s G Cp([d]) by s(k) = ( s(k ) — s(k + 1)) • k 
for all k G [d] and using the convention s(d + 1) := 0. Let there be distributions q,qi,... ,q k G *}3([d]) 
such that Y ^ k =l s(k)Qk = q- Let further U G Mdxd, be a unitary matrix and q G Sp([d]) be defined by 

q(i) := X^fc=i s(&) Y^l=i ^ U ' k ■ Aor ever y k G [d] assume that there exists a p(-\k) G ty(q k ,k). Let finally 
p k G tp([d] fc ) be defined by p k (ii,... ,i k ) := Ke^ <g> ... <g> e ik , U m v k )\ 2 . It holds 

d 

• lb(-| fc ) -Pk\\ > \W~q\\- (32) 

fc =1 

Remark 3. The lemma gains its proper interpretation by letting (ei, pef) = ^2 k= i s(k)\ui k \ 2 be the 
pinching of some state p with spectrum s to the chosen basis. 

The proof of Lemma[2]rests on the validity of the following version of the Cauchy-Binet formula which 
seems to have a certain worth in its own right: 

Lemma 3. Let U G Md be unitary. For every k < d and j G [d] it holds 

k 

K e u ® ...®e ik ,U® k v k )\ 2 =J2\ Uji \ 2 (33) 

(u,...,i fc )e[d]J i=1 


where, as defined in SectionU 3 v k := 'fZreSk s g n ( r )®( r ) e i < 8 > • ■ ■ < 8 > e k . 

Some of our results can only be proven to hold for d = 2, where every representation Vf \ is irreducible. 
These results are summarized below. We start with the definition of two basic building blocks of our 
analysis: 

Definition 8. Let A G £>(C 2 ) and q,p G fp([2]). Then we set 


= -- log(u/M , A^V fM ) 


01 (A,q,p) : 


n 


(34) 






where lim^^oo —/(") = p and — g^ = q. We also set 


e 2 (A,p,q) := - lim -Llog (vf n ,A® 2 - n ,P g( 2 . n) A® 2 - n vf n ), (35) 

n —¥oo Z • Tl 

where by definition of the vectors Vk we have v 2 = i/l/2(ei <8> e 2 — e 2 <8) efi), and finally 

e(p,q,s,A):=- lim ilogtr{P /(n)iA( „ ) ^®"P ff( „ )iA( „ ) At®"} (36) 

n—> oo 77, j i 

where lim^oo ±A (n) = s £ ^([2]). 

Note that @2 does not really depend on p and that the latter parameter is only included into the 
definition in order to be able to deliver a complete description of all quantities within a unified setting. 
The quantities 0, 0i and ©2 are connected via the following theorem: 

Theorem 3. Let lim„_ ) . 00 LfW) = p and lim^^oo 4 g^ = q as well as liniu^oo -A0) = s and A £ Md- 
Then 


©1 (A,P,<l) = H (q) + min D{r\\p 1>A ), 

r£<=-(p,q)} 


® 2 (A,p, q) = ~\ min D(r\\p 2tA ) 

re<$(2,q) 


(37) 


where zz(p,q) := {r £ *)3([2] x [2]) : rq = p, r 2 = q} and pi, A is defined by p\ A (i, j) := |(e^, A'ef}\ 2 and 
P 2 , A {i,j ) := |(e* (g> ej, A^v 2 )\ 2 . In addition to that, 

2 

Q(p,q,s,A)= min V' s(i) ■ QfiA,p, W(Si)). (38) 

W(l»)=9f—T 

l=L 

In the case d = 2 we are able to give a characterization of A u- 

Theorem 4. Let d = 2 and A £ M 2 . The function satisfies, for all p,s £ ^([d]) and <7 £ B(C 2 ), 
the following: 

A A {p,s,a) = —H(s) - s( 2) • logdet AcrA f + s( 1) • D (( p{1) ~ 2 p ) {2) , p(2 ] ( ~* (2) )|[AgA t ), (39) 

where the function D : *p([2]) x M 2 —> R + is given by the convex optimization problem 

D(p, x ) := min ^ P(j) D ( W ( 6 j)\\\ X -A)- ( 40 ) 

W:W(p)=p f—f 

.7=1,2 

Remark 4. Of course this formula demonstrates that A. is continuous in A, p and s - on the region of 
parameters satisfying s{ 2) < p(2). The usefulness of the formula also stems from the fact that it allows an 
explicit and efficient computation of the probability that the state a is detected by a measurement scheme 
which asymptotically detects states with pinching p and spectrum s. 

While it is clear that the function D delivers an efficient way of computing the values of the function A u 
for all unitaries (in fact, for all A £ Md), we have not been able to deliver more insightful reformulations 
of it. We note that interesting connections to matrix scaling (see e.g. the recent work M) are given, 
and another interesting connection is that to information projections in the sense of Wi¬ 
lt seems tempting to look for connections to limc^i D a z for z £ (0,1] but it has already been proven 
that these limits are all equal to D in 12 A . Another possible route would be to look at the limits 
limc-^i -D a r (i_ a ) which we were able to use for r = 1, but we have not been able so far to find any 
relations of these quantities to A a so far. 

At last we exhaust the peculiarities of the case d = 2 to define a set of functions which do, to 
some extent, measure the distance between two states p and a. From the very start, they have offer an 
operational interpretation. In special limiting cases, they deliver either D(p\\a) or lim^^i D a (p\\cr). 

In section eh we provide a proof that they fulfill all the Renyi axioms except the generalized mean value 
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axiom, which is no surprise given that they are well-defined also in the case where they yield the relative 
entropy. 

The fact that they easily deliver lim Q _>.i D a (p\\a) made us step away from attempts to prove that they 
fulfill the data processing inequality. Also, we left any attempts to prove joint convexity to future work. 
We do however prove that the Renyi axioms are fulfilled, except from the generalized mean value axiom. 
First, we need some preliminary notation. Let p, a € <S(C 2 ) satisfy \p, a] ^ 0. Since [p, er] ^ 0 it is clear 
that the set V := {a-p+b-a + c-Id : a, b, c G R} is a two-dimensional real vector space and its intersection 
with the Bloch sphere defines a convex subset of the latter. We may for sake of simplicity assume that 
B is the eigenbasis of p such that 


P = 


5(1 + 2) 

0 


1 ( 1 - 2 ) ) 


(41) 


and 2 G (1/2,1]. We may further assume without loss of generality that the representation of a is such 
that it has only real and positive entries. Both of these assumptions translate to unitary actions which 
depend on p and a. Then, the unitary transformations 


U v := 


cos ip — sin ip 
sin p cos ip 


V G [—7T, 7r] 


(42) 


yield a set of unitary transformations which rotates only the hyperplane defined by p and a. The value 
ip' > 0 at which we get U- v ktU_ , = a-Id + b-p can be used to define the set {EV-t}ie[o,i] which satisfies 
that UqpUq is diagonal and U\aU{ is diagonal in the basis B (the computational basis). If [p,er] = 0 we 
set U t = Id for all t G [0,1]. If d = 1 the convention U t = Idc = 1 applies as well. We are ready for a 
definition: 


Definition 9. Let d G [2], To any p G M ( j satisfying p > 0 and a G <S(C 2 ), set p := trjp}^ 1 • p and let 
Ut := Ut(p,cr) be the set of unitary transformations which arise from p and a as described above. This 
defines a set {-Rt}te[o,i] of relative entropy like functionals via 

Rt{p\W) ■= Ayt (phich(£/tp[/j), spec(p), tr{p -1 } • a). (43) 

Remark 5. Theorem Q] ensures that the definition does not only lead to trivial concatenations of ro¬ 
tations followed by unitary transformations, since Rq = D 1 and R\ = D are different functions. The 
normalization factor tr{p -1 } in front of a enables one to prove the order axiom. 

The structure of the functions we defined so far delivers operationally meaningful quantities right 
from the start, as they describe the asymptotic scaling of the probability that certain tests yield specific 
outcomes given that a system is in state er®" and n is large. 


IV Proofs 

We now give the proofs of our theorems, in order of appearance. 

Proof of Theorem [3 The proof of statement 2 is implicit in |20l and what is left to do is giving the proof 
of statement 1. Any of the representations Vf t \ is irreducible if f(i) = A T (i) for some permutation r G Sd- 
This can be seen as follows: Denote the set of all tableau T of shape A by T*. Then statement 1 can be 
seen to hold true as follows: Remember that each V/ is invariant under B, so that for each T G T>, we 
have EtV G Vf whenever v G Vf, where 

Et := E E sgn(u)B(i> o t). (44) 

V(zlCt t£Rt 

is the Young symmetrizer corresponding to the tableau T, Rt is the set of permutations which permute 
only the elements in each row of T amongst each other and Ct permutes only the elements in the columns 
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of T. By [221 Chapter 5.5] (replace the object T n C d := (C d )®" there with Vf) it holds that for every 
fixed T £ T\ the dimension of the vector space spanjEri; : v £ Vf} gives the multiplicity of F\ within 
Vf. Now let for sake of simplicity / = A. It is clear that for the vector v = and T be ’the’ 

standard tableau with numbers 1,... ,n filled in starting from left to right in the first row, then carrying 
on from left to right in the second row, and so on. 

Then v := E T v ^ 0. Thus dim V\.x > 0. Now take any other product vector w = 0" =1 e Xi w here x n £ Tf. 
There is at least one column (say the first) having at least two equal entries (for example it holds that 
both x\ = 1 and x\^+\ = 1). This statement is valid as well for every w T := B (r)w whenever r £ Rt, 
only the position of the specific column changes. Take the permutation 7r = (l,Ai + 1) £ S n which 
interchanges the elements X\ and x\ 1+ \ for every x n £ [d] ra . It holds B(7 t)u; = w. On the other hand 
7r £ Ct and sgn(7r) = — 1. For each r £ Rt , let i r T £ Ct be a corresponding permutation that satisfies 
B(7r T )w T = w T and sgn(7r r ) = —1. It follows 

Etw = E E sgn(u)B(i))ui T (45) 

= E E sgn(u)B(u) • B(7r r )u; T (46) 

T(z.Rt V(z.Ct 

— EE sgn(u o r)B(u) • B(7r r )w r (47) 

= - E E sgn(u)B(u)w T (48) 

v(z.Ct 

= — Etw , (49) 

so that dim = 1 follows. The argument is independent under a transformation / i —> f o r whenever 
r £ Sd, so that all the representations Vf t \ for which f(i) = A r (,) for all i € [d] holds true for some t £ Sd 
are irreducible. 

Now we connect our first observation to a trick that we shall use more often in what follows: 

Let V be an irreducible subspace of the symmetric group. Then for every 0 ^ v £ V we have A v := 
reSn ^iB(r)|u)(u|B( T_1 ) = c • Py for some c = c(v) > 0 and the orthogonal projection P v onto V. 
This is seen as follows: note first that A v ^ 0 whenever v ^ 0. Furthermore each A v is invariant under 
permutations. By Schur’s Lemma (see e.g. [22] Chapter 2.3]) it follows that A v = c(v ) • Py. 

Moreover, by taking the trace we see that ||u ||2 = c • tr {Py}. This implies that for every a £ C d and 

v £ V we have 


tr {P v a® n j 


tr {Py} 

Ml 

dim(C) 

TIT 


tr{|u)(u|a® n 

-(v,a® n v). 


} 


(50) 

(51) 


Thus all that is left to do in this case is to construct one vector v within Vf x and calculate its norm as 
well as (v, a® n v). 

This task again is straightforward since we may just use the standard tableau T that we defined already 
and the corresponding Young symmetrizer Et- Applying this symmetrizer to the vector yields 

(without loss of generality / = /■*•): 


v := E t ®? =1 ef /(i) 

= \Rt\ J2 sgn(r)B(r)®t ie f /(i) 
tGCt 

= |7?t|B(t / ) <8)f =1 ( sgn(r)B(r)ei 
VreSi 


) 0(Ai —Ai+i) 


(52) 

(53) 

(54) 
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where t' is a suitably defined permutation. For any k £ [d] we may now define v k £ (C d )® fe by 


v k ■= sgn(r)B(r)ei <g) ... <g) e k , 

TeS[fc] 


( 55 ) 


a shorthand that allows us to write 

v = |i?T|B(r') ® d =1 wf (Ai_A<+l) . (56) 

In order to get a lower bound on the norm of v we first note that each n, satisfies \\vi\\\ = il, so that 


= i^t| 2 II(* ! ) 


Ai-Ai. 


Another important ingredient is the equality 


E E sgn(r) sgn(cr)B(r)B(u) = k\ ^ sgn(r)B(r) 

reSfc v&Sk T£S k 


which lets us conclude that 


We are finally able to compute 


(v k ,a® k v k ) = k\ • det(c7i :fci i :fc ). 




INK 




i=1 
d 


^E|i? T | 2 n(*!det(a 1:! ;, 1:i )) Ai - Ai+1 

11^112 i=1 

^SrVl 2 (ri(i!)—) (n 




dim(V’) 


v?=i 
f d 


\3 =1 
d 


\ rt \ 2 Ai+i n aet (-ww)) 




Ki=l 


U = 1 


= dim(V) <let Sdr; ;; i ■ i )) Aj Aj+1 

3= 1 

= pl( ? x)2 rl ' ff ( A )2"'^-i ! = 1 ^ Ai_Ai+1 ^ logdet ^' Tl:i ’ 1: d 

= pl(?i)2 n ' ff ( A ) + ^- 1 ( Ai - Ai + 1 ) losdet ( CTl ^ 1: d 
_ pl(n)2 n ' ,s '( p ) + ^-'*= l( ' Pi_Pi + 1 ) lo s det ( <T i^.i:i)+e(n) 

= pl(n)2 n ' (£l (pll d )+dn)) 


(57) 


(58) 


(59) 


(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


where d is a nonnegative matrix which is simultaneously diagonal with p and is defined via its diagonal 
entries an = det^i^m) and lim^-yoo e(n) = 0. In [1] (see Theorem 2 with the respective parameter r of 
the theorem set to r = — 1 and equation (24) there) it has been proven that D(p\\a) = lim a _n D a (p\\a), 
so that ultimately we have 

— lim — logtr{P/ act®"} = lim D a (p\\a) (69) 

as desired. □ 
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We now start to investigate the scalar products 


tr{A® n P ft \A*® n P f i j \i}. (70) 

It is generally clear that A = A' has to hold, so that the questions we pose get reduced to the evaluation 
of quantities of the form tr{A® n P g A 1l ® n P\Pf}. 

Proof of Theorem [H We now work with an arbitrary d £ N. We write 

tr {A®-P g A^P KX } = ^|^tr{|t,><<;|P g } (71) 

as before, and again v takes the form 

d / \ ®(Ai —Ai+i) 

v = |-Rr|0) [ ^2 sg n ('r)B(r)e 1 ® ... e* J . (72) 

i=1 VrGSi / 


Again we set Vk := X^reS fc sgn(r)B(r)ei <8 >... <8> e*,. The asymptotic behaviour of the function (/, A) i-» 
— logtr{P/^} is known to equal that of A >->• Alogtr{PA} for all / and A satisfying Kf t \ > 0, so that 
what is left to do is the following: We have to calculate 


d 

tr {A® n P g A^® n 0 |n l )(u i | 0(Ai - Ai + l) } 
2=1 


d 

J2 J|tr{A t ®( Ai - Ai + 1 )p ffi A t ®( Ai - Ai + 1 ^ l )(w i |®( Ai - Ai + 1 )}, 

gi+...+gd=g i =1 

(73) 


and again we now have to dive into calculating, for every to £ N and h £ T m and t £ T m >, as well as for 
every k, quantities like tr{A® k ' rn PhA^® k ' rn \vk){vk\® m }- This task needs some additional notation. Let 
to = in' ■ k for some natural numbers m' and k. Then, we set 



It then holds that 


if E(u,...,u)e[rf]^(*i>--->4) = h{j) V j £ [d], 
else 


(74) 


l Th > J2 (75) 

t£T m ,([d] fc ) 

Moreover, those types t which do not occur on the right hand side of inequality d75l) but only on the left 
are exactly those which have t((ii,..., ik)) > 0 for some choice (*i,..., ik)- Such types however satisfy 

(v k ,A® k e il ®8 e ifc ) = 0 (76) 


by symmetry of Vk- This justifies (actually it does so only in the second row of below chain of estimates 
so one has to read from there both back- and forwards) that we write 


tr{A® m - k P h A^® rn - k \v k )(v k \® m '} = ^F/Tr{A® fc PiA t ®^> fc )(v fe |® m '} 

t 

= J2 F t h II (vk,A® k e il ® ...®e ik ) t{i "-' ik) 

1 (u,■■■,*&) 

_ ^ ph2 m '-T2(i 1 ,...,i k ) 7^7*(*i>---dfe)log(|(uf,,eiii8i...ei fc >| 2 ) 
t 

< pl(n)2 m, ' maXt t(n,...,u)log(K^,e i i0...e ifc )| 2 )) 


(77) 

(78) 

(79) 

(80) 
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Upon normalization, the definition of Fj 1 translates into the set of probability distributions on [d] k which 
we defined in the introduction: The set 




p e W] fc ) 


... ,i k )) > 1 => p((ii, ■ ■ ■,h )) = 0 
P{[d\i) = k ■ q(i) V i € [d\ 


(81) 


It is this set that determines the asymptotic behaviour we are after: namely, for lim^-nx, A h= q it 
holds that 


lim —log tr {P h (k-n) \v k )(v k |®"} 

n-¥ oo K ' Tl 


1 

— min 

k pe<#(q,k) 


D{p\\pk) 


(82) 


with p k € fP([d] fe ) being defined viap fc (ii,... ,i k ) := Ke^®.. .®e ifc , A® k v k )\ 2 for all (i lr ...,i k ) e [d] fe . For 
sequence of partitions g [ n ' + ... + g^ = g^ (for each g ["' l we have n k := Y^=i 9 (j) = k • ( X k — Afc + i) 
and lim„_ > , 00 n k /n = s(k) for all k = 1 ,..., d) with respective limiting distributions q, q ±,..., qd we get 


lim -logtr{P ( n)P A( n) A (n)} = lim - log tr{P A(n ) }tr i P ( n> |u fc )(u fe |®( Afe Afc+1 ) l (83) 

n—¥ oo 77 ^ n—>• oo 77 

Kk=l ) 

= H(s) + y: lim - logtr {p („)|u fc )(u fe | (8l ( Afc _Afc+1 ) 1 (84) 

z ' n-±oo Tl y k 

k—1 K J 

= H{s) ~ U w* », D < 85 > 

fc—1 


Naturally, this leads us to the formula 


lim - log tr{P ( n) P A( n) A( „)} 

n—>■ oo 77 y 


lim — log tr{P A (n) }tr 

n—>oo 77 




/c=l 


= P( S ) 


d 


min > s(fc) min 

peq3(w((5 fc ),fc) 


• -D(plbfc) 


( 86 ) 

(87) 


which is valid for all asymptotic shapes s = lim„_>. 00 -i-A^ and where s(k) := (s(fc) — s(k + 1)) • k and 
s(k + 1) := 0. Note that 


d d d 

s ( fc )= ■k-J2 s ( k+1 )' k 

k= 1 k=l k=l 

d d 

= s • k - H s ( fc ) • ( fc - !) 

fc=l k—2 

d 

= s{i)+J2 s ^ 

fe =2 


Let now A= U G Md be a unitary matrix and (q k ) k=1 and (p(-\k))^. =1 be such that 

d i 

6 (q,s,U) = H(s) - ^2 s(k)-D(p(-\k)\\p k ). 

k =1 


( 88 ) 

(89) 

(90) 

(91) 

(92) 
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By Pinsker’s inequality and convexity of x H > x 2 we get 


&{q,s,U) < H(s ) - 2^Ts{k)l\\p{-\k) - p k \\ 2 

k= 1 

( d \ 2 

5Z J ( fc ) Vfc ll-P(‘l fc ) -Pk\\j ■ 


( 93 ) 


(94) 


It is trivially true that — 4 < — 4 for all 1 < fc < d, so that by monotonicity of the square root we get 


®{q,s,U) < H(s ) - 2'}Ts(k)4~\\p(-\k)-p k \\ 


k= 1 


< H(s) - \ (^2s{k)\\p{-\k)-p k \\ 


(95) 


(96) 


v.k=i 


The distributions qi,... ,qd and p(-|l),... ,p{-\d) fulfill the assumptions of Lemma [2J so that by convexity 
of || • || we get 


0(q, s, U ) <H(s)~ |||q - <?|| 2 , 


(97) 


where q(i) := Ylt -i ®(^) i f° r ad * €= [d]. Of course then, Q(q,s,A) attains its maximum H(s ) 
when q = q. This is the only maximum, since the function q K > O(g, s, [/) is convex: Let A £ [0,1] and 
set A' := 1 — A. Then for arbitrary q , g' £ f£([d]) satisfying s ^ q and s ^ g' we have 


' 1 

0(Ag + X'q', s, U) < min min ) s(k) — min D(p\\pk) 

' ~ ff(s)=,r( s >«' ^ k P G(\w+\'W')(6 k ) " ' 

/C—1 

^' 1 

< min min N s(fc) — min min D(Xp + X'p'llpk) 

W(s)=qW’(S)=q' k p£W(,6 k )p'£W(5 k ) 

d 1 

<min min Vs(fc)- min min. (A-D(p||p fe ) + A'D^'llp*,) 

W(s)=q W'(s)=q' ““J /C p£W(<5 fc ) p'GW(<5 fe ) 

^ ^ 1 

= min min Vs^-A min D(p||p fc ) + A' min ; fl(p'|W 


k=1 


p'ew(8 k ) 


d d N 

= min min ( ^r^A min -D(p||pfc) + ^y^A' min -D(//||pfc) 

wcsw„w/rsw n / \ ^ k p & w(5 k ) ^ k P 'GW(S k ) J 

\k—l k— 1 / 


W(a)=q W' (s)=q' 
d 


= X min Vf min D(p\\p k ) + X' min V] ^ min L>(p'||p fc ) 

VF(s)= 9 nCW(A. ) M/'f.S \=n> ^ k n 'GW(X,.\ 


k= 1 


pew(s k ) 


Vy / (s)=g / 


fe=l 


p'ew(8 k ) 


= XQ{p,s,U) + X'G(p',s,U). 


(98) 

(99) 
( 100 ) 
( 101 ) 
( 102 ) 

(103) 

(104) 


□ 


We now give the proofs of our two additional Lemmata. The asymptotic estimates for 
tT{A® n Pf'xA'® 11 Pfi t \/}, although they give a rather cumbersome impression, naturally introduced the 
distributions p k £ fp([d] fc ) defined by pi(ii, ■ ■ ., ii) := Ke^ ® d k , U® k Vk)\ 2 - These again make it interest¬ 
ing to look at lower bounds on the exponent in terms of norms, which turn out to deliver a satisfying 
intuition once Lemma [3] holds. We will now prove this Lemma. 


15 



Proof of Lemma\^ Let us assume that j = 1 holds. It will become evident from our proof that this 
is without loss of generality. We will have to consider appropriate submatrices of U , that are defined 
entry wise as follows. 

U[(il)Ll,(jl)Ll}mn--=U imjn . (105) 

For k a natural number, the letter k denotes the string (1,..., k). For in £ [fc], k\m denotes the string 
(1,..., to — 1, to + 1, k). We let [d](2, k) be the set of all strings of length k with elements taken from 
{1,..., d} without repetition. Using these matrices will allow us to employ first Laplace’s formula, then 
the Cauchy-Binet formula followed by the Sherman-Morisson formula. Together with the fact that U is 
a unitary matrix, this will lead to the desired result. 

Now, we will write above sum as a sum over determinants. This will allow us to apply the Cauchy-Binet 
Formula - but first we have to rewrite our form slightly in order to see the determinants. 


Y\( Ue n ® ■ ■ • Ue ik ,v k )\ 2 = ^|(Ue 4l (8).. .Ue ilc , Y sgn(r)e T( i) <g>... <g> e r(fc) | 2 (106) 

[d\$ [d]J v fc reSfc 

= l E 1 E s g n ( r KiUi) ® ® u i kT {k) | 2 (107) 

[d]J reSk 

= ^|de t ([t/,(* i ,...,i fc ), k ])| 2 . (108) 

Mf 

It is clear that the terms in above sum are invariant under permutations. The function (i\,... ,ik ) Ht 
KC/eij <S> ■ ■ ■ Uei k ,Vk)\ 2 is also designed such that \{Uei x ® ... Uei k ,Vk)\ 2 = 0 whenever i m = i n for some 
to 7 ^ n. This implies that it suffices to consider those terms where 1 stands in the first place and 
(1, * 2 , - - -, ik) form an index set (meaning that [i 2 ,... ,iff) £ [d](2, k)). We thus get the formula 

Y\(Ue il ®...Ue ik ,v k }\ 2 = £ |det([U, (1, i 2 ■ ■ ■, ik), k]) | 2 • (109) 

[d] f (i 2 ,...,ife)e[d](2,fc) 


Since only the columns 1 to k enter our calculations, let us consider U as a d x k matrix from now on, 
with the transposed matrix U T being a kx d matrix. U denotes the matrix having the complex conjugate 
entries of U. We will now apply Laplace’s formula (twice), followed by the Cauchy-Binet formula [Tdl 
Chapter 0.8.5] (set r = k — 1 in the book): 

D<c e ix <g).. .Uei k ,Vk)\ 2 = E |det([t/, (1, i 2 ..., ife), k])| 2 (110) 

[d]f (*2,...,ife)e[d](2,fc) 

k 

= Y E (-l) m+n ui m uUfdet([/7, (i 2 ... ,ifc),k\TO])det([i7, (i 2 • • • ,4), k\n]) (111) 

(<2,... l <fc)GH(2 l fc) m,n= 1 
k 

= Y E (-l) m+ "«imUtedet([U T ,k\TO, (* 2 ---, 4 )])det([L r , (i 2 ■ • • ,4),k\n]) (112) 

(* 2 ,-.,ifc)G[d](2,fe) m,n=1 
k 

= Y (-l) m+n Mi mUm Y det([C/ T ,k\?n, (i 2 ...,ifc)])det([U, (* 2 ...,ife),k\n]) (113) 

m,n=l (i 2 ,...,ifc)e[d]( 2 ,fc) 

k 

= Y, (-l) m+n Ui m Mmdet([U T ,k\TO,d\l] • [U. d\l, k\n]). (114) 

m,n= 1 

It looks tempting to re-apply the Laplace formula here, but the determinants are now being calculated 
on products of non-square matrices so that we have to find a different means of dealing with the above 
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sum. Let us calculate above determinants. For a fixed pair (to, n) £ [k] x [k] we have to calculate the 
determinant of the (k — 1) x (k — 1) matrix C(m,n) defined by 

C(m, n) := [U T , k\m, d\l] • [17, d\l, k\n]). (115) 

The entry of this matrix that corresponds to * £ [/c]\{to} and l £ [fc]\{ra} is given by 

d d 

55 [M d\l, k\m]), r ([?7, k \n, d\l]) ri = 55 (116) 

r —2 r =2 

= S(i,l) - uuuii, (117) 

and it is exactly here that we use the fact that U is a unitary matrix. Since all entries belonging to 
the TOth row index and the nth column index are removed from C(m, n), it is evident that det C(m,n) 
equals the (to,?z) minor of the k x k matrix ( S(i,l) — uuuu)u, which can equivalently be written as 
M := 1 — \u)(u\ £ Mk where u = JT=i u ii e i- The determinant of M is calculated as 

k 

det(M) = 1 - 55 |u H | 2 (118) 

i=1 

via [9J Lemma 1.1]. This makes it useful to again apply Laplace’s formula (twice, again), where M mn are 
the entries of M: 


k 

1 - 55 M" = det(M) (119) 

i= 1 

k 

= 1 (-l) m+n MnmdetC{n,m) (120) 

m,n=1 
k 

= i 55 (-1 ) m+n (S{7n, n) - -»i m ui ra )detC , (n,, to). (121) 

ra,n—1 

It follows again from El Lemma 1.1] that for every to £ [fc] we have detC(m, m) = 1 — X5=i l^iiP + l^iml 2 , 
so that 


k 

i - 55 m 2 = 
2—1 


1 + *:55m 2 - 

k 

-j; 55 (-l) m+n Mi m ui ri detC'(TO,n) 

m,n=l 

(122) 

= 

(i-ElM 2 ] 

1 + 1 55 m 2 - 

T, 55 KM ® • •• ®Ue ik ,v k ) | 2 , 

(123) 


V »=i / 

2 — 1 

H{ 



which is equivalent to saying that 

k 

55 m 2 = 55 Kt/e^ <g>... ® Ue ik ,v k )\ 2 (124) 

i=i [d]f 

so that we have proven the desired formula. □ 

We now come to the proof of the estimate which originally motivated us to study the determinant 
equation which is the content of Lemma [3] 

Proof of Lemma\^ Let the preliminaries of the Lemma be fulfilled: We have that s £ fpl([d]) and 
define s £ *}3([d]) by s(k) = (s(k) — s(k + 1)) • k. Let there be distributions q, q ±,..., q k £ fp([d]) such 
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that Ylt=i s(k)hk = <?■ Let further U £ Mdxd, be a unitary matrix and q £ *P([d]) be defined by 
q(i) := X)fc=i s(fc) Y2i= 1 • F° r every fc e [d] assume that there exists a p(-|fc) £ fp (qk,k). Let finally 

Pk £ < P([d] fe ) be defined by p k (h, ■■■ ,ik) ■= Keq ® ® e ik , U® k v k )\ 2 . Let [d]% R := U 9 e<p([d])^P(9> *0- 

These are those sequences of length n with elements taken from [d] that have no single element occurring 
twice (NR means “no repetitions”). It holds 


d d 


^2 s(k)\\p(-\k) - p k \\ = ^2§(k) ■ 
k =1 k= 1 

E b(*i--- 

• ,u|fe) -Pk(h, ■ ■ 

* 5 | 

(125) 

d 

= H k ) ■ 

k= 1 

E 

(ii...tis)e[d]5vR 

• ,u|fc) -Pk(h, ■ 

• • 5 | • 

(126) 


Since no repetitions are allowed in above sum we can be sure that, for every k £ [d], the respective sum 
over [d\ RrR can be split up into sums over subsets [d]£ as follows: 


E \p(ii,...,i k \k) 

(i 1 ...i k )&[d]% R 


d 

Pk(h, ■ ■ ■ dk)\ = E 

i= 1 


l E 


\p(ii, ■ ■ ■,ik\k ) 


Pk(iu---,ik)\, (127) 


where the fact that [d]* (~1 [d]| ^ 0 has been taken care of by the factor p We can use this to reformulate 
the above sum as 

d d^ 

j:mM-\k)-P k \\=j:m V E \p(ii,---,ik\k)~pk(h,---,ik)\ (128) 

fc=l fe=l (ii,...,i fc )£[d] fc 

d d 

= '52( s ( k ) - S (k + 1 ))'Y1 E \p(k,---,ik\k) ~Pk(ii,---,ik)\ (129) 

k=1 i=1 (ii,...,i fc )e[d]f 

d d 

>^2( s (k)~s(k+i))j2 1 E (p(h,..., i k \k) - p k (ii, ■ ■ ■, ik )) | (130) 


fc=1 i=1 

d d k 

= E w fc ) - < k +!)) E i fc ) - E m 2 i ( 131 ) 

k= 1 2=1 Z=1 

d d k | |2 

= E ( s ( fc ) “ s ( fc + x )) E ■ 9fc(*) - k ■ E (132) 

fc=l 2=1 Z = 1 

>DX>)-«M-X>)£bA ass) 

2=1 k=l k=l 1=1 

d 

= El 9( l)^(*)l (134) 

2=1 

= 119-911- (135) 


□ 


V Proofs for d = 2 

Proof of Theorem [^} In order to go further with our investigation of the asymptotic behaviour of 
tr{P/^cr® n } for arbitrary / and A we unfortunately have to live with the restriction d = 2. 

In our case there is no difference: While it seems to be a rather involved task to obtain explicit formulas 
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for the case / ^ A whenever d > 3 we are well able to so when d = 2. The reason for this is that in this 
case we always have Kf.\ £ {0,1}, so that each Vf t \ is irreducible ! This can be seen as follows: 
According to [22J Chapter 5.5] (replace the object T n C d := (C d )®" there with Vf). the multiplicity of F\ 
within Vf t \ is given by span{£Yv : v £ Vf} where T £ T A is any standard tableaux of shape A. 

Let again T be “the” standard tableau with entries Ty = (j — 2) • i + (j — 1) • (Ai + i) and w = 0" =1 e Xj 
for some x n £ Tf. Then 


Etv = E E sgn(u)B(u) -B(r)u. 

t£Rt v£Ct 


(136) 


Define At := J2 v eC T s g n ( v )®( w ) and Bt := ^2 t€Rt B(t). For every r £ Rt , set w T := B(r)v. For every 
r £ Rt we see that Atw t = 0 holds if and only if rx n = (x T -in\,... ,x T -i („)) satisfies (rx n )i = ( rx n )j 
for some pair (i. j) where 1 < i < X 2 and Ai + 1 < j < n. It follows that Atw t = c(v, t)v for some set of 
non-negative numbers {c(r), r)} Te 5 n and a vector v to be calculated more explicitly later. Thus 


Etw = 



(137) 


proving that the multiplicity of F\ in Vf t \ is at most one. We proceed with the calculation of Etw. 
Again, take the Young symmetrizer Et x ■ Then for some constant c'(/, A) we get 

v := Er x eU ef«-» 

= c'{/, A)B(t')(T„ 2 )®(*.-M g, 

where Vf-\ 2 is defined only for those pairs (/, A) for which both /(1) > A 2 and /(2) > A 2 holds. In that 
case / — A 2 defines a new type g := f — A 2 on {1, 2} Al-A2 so that it generally makes sense to define for 
an arbitrary g £ T m : 

v g ■= Er E ( 14 °) 

x m eTg 

The asymptotic scaling of tr {Pf \a^ n } is then conveniently calculated by starting with 

-Eogtr{P /iA o-® n } = -Eog (^^c(f,X) 2 (^{v 2 ,a 02 V 2 ) X2 {vf-x 2 ,a 0{f - Xll) Vf-x 2 )j (141) 

= (logdim(Y A ) + log(i(u2,cr® 2 ^2) A2 ) + log((u / _ A2 ,cr® (/_A2) u / _A 2 ))j , (142) 

then calculating the limit of the three terms in the sum separately yields the desired result - but only 
if the limiting behaviour of the last one of them is known. We thus start with that part. Under the 
assumption that (4-/") n e n converges, we call the limiting object p := lim,!-^ A/". For any given state 
<t on C d with matrix representation a = JY . <Tj j\e-i){ej\, it is then of interest to describe the limit 

lim —log ((uf™,cr® n u/n)). (143) 

n—>-oo n 

This can also be cast in to the form of the subspaces V/ iA by setting A = (n, 0,0,..., 0). In order to have 
a more streamlined notation, we will drop the superscript n in f n for now, then we can upper bound the 


(138) 

(139) 
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limit as follows: 


1 n 

{ V f,a® n v f ) =— nw™ 

1 / 1 x n ,y n eT f i—l 

(144) 

71 

(for some fixed y n £T f ) = ^ Jpcr ) XitVi 

(145) 

x n eTf i—l 


= £ 

(146) 

x n €Tf i—l 


2 

_ ^ JJ c N{i,j\x n ,y n ) 

(147) 

x n £Tf i=l 



where we have set := Now for every pair x n , y n £ Tf it is clear that the numbers N(i,j\x n , y n ) := 
\{k : Xk = i, Uk = j}| satisfy N(i,j) = N(j,i) for all i,j £ {1,2}, this being another peculiarity of the 
case d = 2. Note that this implies 


(v f ,cr v f ) = 


n 


„N(i,j\3 


x i. -i 1 A. "V 

GT f i=1 


(148) 


Obviously we need some additional structure. This comes into play by decomposing the set Tf according 
to 


Tf = U G LA X T g2 

9i:92 


(149) 


where each T 9i C [2]^W, and the above union is over disjoint sets. The numbers £ {0,1} are 

defined in analogy to the Kronecker coefficients of the symmetric group, precisely speaking we set 


f 1, Ej=i&0) = /(*) V i £ [d\, 
\ 0, else 


(150) 


The nice thing about this decomposition is that there are only polynomially many (in n ) different choices 
( 51 , 52 ) - more accurately, the number of such choices can be given a the loose upper bound pl(n) := (2n) 4 . 
This allows for a reasoning along the lines of the ’method of types’: 

<»/.»*"»/> = V C'n \T„\ ■ eg" 1 (151) 

91,92 i,j=l 

d 

<pl(n) max Yi \T gj \ ■ c 9 ^. (152) 

Si,S 2 :Gj 1 , 92 > 0 ij . =1 

As a consequence of an almost identical calculation it follows that 

2 

(vf,a® n v f )> max \\ (153) 

gi ,92 :G 91 )92 >0 i j = \ 

This demonstrates that the following holds: If limn^oo —= p, then 


lim 

n —>00 


n lo S( V T 


.*® n vr 


) = lim 

n —>00 


■log 


max 

91 ^ 92 'Gg-^ t g 2 



(154) 
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holds whenever the two limits exist as well. We now translate our statements to a different regime by 
noting that f{i)~ 2 \T gi \ < 2 nH ^ i \ where Qj := We then estimate 


{v f ,a® n Vf) < max 

5oi-■ ■ i9d'G 9 Q t gi >0 

= max 


(155) 

(156) 


If we now plug in the limiting behaviour (■^■/ TI ) n eN —> P and translate definition 11501 to probability 
distributions by dividing through n, we end up with 


E G L»n E IK™ 

(157) 

51.32 3=1 yxf^GTg. i —0 / 


< max 2 _n ^i=i p(j) £> ( w ( 5 i)ll c jl) 

(158) 

W:W(p)=p 

— 2 ~n min w , W (p)=p I2j= 1 p{ 3 ) D ( w ( s j)\\ c -j\) 

(159) 


and the symbol W stands for the matrix (®ij)y =1 with nonnegative entries and u;(l|z) + w(2\i) = 1 
for i = 1,2, and W(p) = j=i w (i\j)p(j)^i can be seen as application of the matrix W to the vector 
p = where 5i(j) := S(i,j ) are the usual Dirac distributions on [2]. 

The calculation of a corresponding lower bound can be established with an almost identical reasoning, 
so that we obtain 


J™, n 1 °g<u/",CT® n v / ») = - min ^p(j)D(W(5 3 )\\\a. 3 \). 

n—r oo w: W (p)=p — 

3 =1 


(160) 


We collect what we found so far in the following formula: For spec p = (p\ , 9,2) with p 1 > P2 and pinching 
pinch p = {y\, 1/2) we have 


$0lk) = S(p) + p 2 logdet(a) + {p x - p 2 7 I 7 ^|)l|cr)- 


(161) 

□ 


We now turn our attention to the scalar products for the special case d = 2, which allows for some 
stronger results. 


Proof of Theorem 0 Using the same tricks as in the proof of Theorem [5] or Theorem [T) we can write 


tr{Pf,\A® n Pf> A^® n } = 


dim Vf t \ 


tr{A® n \v)(v\A® n P f , } y} ■ S{ A, A') 


(162) 


where v is as defined below and with respect to a standard Young tableaux T that we write e.g. for 
n = 15 as 


1 

3 

5 

7 

9 

11 

13 

15 

2 

4 

6 

8 

10 

12 



such that the role of the anti-symmetrizer Bt = Y1v£R t sgn(u)B(u) is to anti-symmetrize on the first 
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2 • \ 2 blocks. We can then write (setting g(i) := f(i ) — A 2 ): 


1 n, ^2 

^ /2 VW7\ T 1 

(163) 

2 A 2 /2 mn E E sgn(u)B(u • r) ®? =1 ef /( } 

Z VHsI T£R T veC T 

(164) 

oa 2 /2 ( E sgn(v)M{v) 0(ei ® e 2 )^j 0 v„ 

\v£Ct i= 1 / 

(165) 

vf A2 0u ff , 

(166) 


and what remains is to calculate the quantities ||?;|| 2 and ( v,P'j,v }. The former evaluates to ||w||| = 
2 Aa • \T g \. We calculate the latter by exploiting the specific product structure of v that we developed 
above. Going into details, we see that 

(v,A® n P r A*® n v) = J2 G f g ' ug2 ti-{A® 2 - X2 P gi Ai® 2 - x *\vf X2 ){vf X2 \} ■ tr{A^ Xl ~ x ^ P g2 

91,92 

(167) 

holds. It is now clear that we have to calculate, for every even natural number n and type g £ T n , 
the quantity ti{P g A® n vf n/2 Al® n } and for arbitrary n£ N the number tr{A® n P g Al® n \v g ){v g \}. This is 
done in the following. We define a function p 2 ,A ■ [2] x [2] —» R+ via P 2 ,A(i,j) ■= |(e* ® e 2l A® 2 v 2 )\ 2 (if 
A is a unitary matrix this is an element of fp([2] x [2])). Again we look, for every two types g £ T„([2]) 
and h £ T n / 2 ([2] x [2]), at the numbers that we defined in d74l) . These enable us to write 


tv{A® n P g A^® n vf n/2 } = Y {x n ,A'® n vf n/2 A® n x n ) (168) 

x™eT g 

n/2 

= e ni< e -® e - + ^ t ^)i 2 ( i69 ) 

x n £T g i= 1 

= Y F h\Tk\ n \{ ei ®e h A^® 2 v 2 )\ 2 - h ^ (170) 

h i,j =1 

< pl(n)max2"l (H( ' i)+i: ^=^ (i ’ J ' )logP2 ^ (i ’ J ' )) . (171) 

h 

The estimate can be carried out in the other direction as well such that we get, for every sequence g^ 
such that lim^oo 1g = q £ *P([d]) holds, the asymptotic relation 

lim — logtr{A® n P g A' i ' 0n |w®"^ 2 )(i;f ln ^ 2 |} = — 1 min D(r\\p 2 a), (172) 

n -¥00 n r£<$(2,q) 

and this does obviously imply that 

0 2 (A q) = ~\ min D{t\\p 2 ,a)- (173) 

r£<V(2,q) 


The other asymptotic quantity that needs to be calculated is still left open. Here, we proceed as follows: 
Define 


r<g'-yg ._ / 1> Si=i 9j{ *) — d{j)> i — 1,2 
91,92 ' \ 0, else 


(174) 
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where each gj £ then it holds that withpi iA : [ 2 ] x [ 2 ] —> K + viapi iA (*, j) '■= |(e*, Vej )| 2 (pi,a('j*) 

is a an element of t)l([ 2 ]) if A is unitary) that 


( V : 


3 ’ ' 


% P 

- 1 n 


\) <pl(n)2 n ^ >\(e x n,A^ n v g )\ 2 

= pl(n)2"^')-*(s)| £ (e x „,At®- e?; n )| 2 


S”GTg 


= pl(n)2 nff( 3' ) - ff( 3 ) • 


E g s^Mt 5 iHt 52 |. 


Sl .32 


*,j=l 


(175) 

(176) 

(177) 


< pl(n) max 2 n ^ H ^ 9 '^~ H ^ +2 ^= l9 '^ H( - 9i ' ,+ ^’i= l9 '^ 9i ^ los ^ ei ’ Afe:i ^\ (178) 

_ Sl+S2=S 

and since an equivalent lower bound can be established as well we obtain that for two sequences (</^) n eN 
and (<?( n ))neN with respective normalized limits q and q we will have 


lim. K(«) ’ P'gM v g w) = -H(q)+ 


(179) 


+ max S2q(i) \ -\ogq{i) ~'S2w{j\i)logw{j\i)+ 'S2w{j\i)logpi :A (.iJ) 

wmv=y ^ 

(180) 

2 

= -#(<?) + max E 9 (*)^(j|i) (-logg(i) - log«j(j|i) + logpi iA (*,j)) (181) 

W:W(q)—q r — 


*,i=i 

= --ff(g) - min_D(r||pi iA ), 


rGH 


(182) 


where S := {p £ *p([2] x [2]) : pi = q, p 2 = q}. Therefore 

@!{A,p,q) = H{q) + min D(r\\pi, A ) (183) 

rGS 

We have thus identified the building blocks of (v, P'pv) and can now write 

0(p, q, s, A) = w mm ? (s(l)0i(A, (fgEfjg, WW) + S(2)Q a (A, W(S 2 ))) (184) 


□ 


VI Axioms 

We give a short overview over elementary properties that make the Rt candidates for relative entropies. 
We note again that various other possibilities exist to define one-parameter families of unitary transfor¬ 
mations - e.g. via utilization of the geodesic (see [4]) t >->■ 7 Pi<T (f) where 7 p , a (t) := p 1 / 2 (p _ 1 / 2 crp _1 / 2 ) p 1 / 2 
which (upon normalization) draws a path between p and cr that (like the definition that we use here) 
enables one to uniquely define a pinching of p to the eigenbasis of 7 (f) whenever [p, cr] = 0 holds. Since 
the geodesic curve obeys ^upW ,UaW (t) = U”fp,cr(t)W for all states p, a and t £ [ 0 , 1 ] and unitary 
transformations U this definition leads to another quantity, call it R t , which is unitarily invariant just 
like R t was. We leave further investigations of these connections to future work and look at some 
properties of the Rt family: 

Continuity. Continuity follows directly from the fact that both functions can be rewritten as a 
convex optimization problem where both the function to be optimized and the convex set that it is being 
optimized over depend continuously on p and cr. 
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Unitary invariance. The choice of the basis B in which the Pf.x are defined is just such that 
it changes as B a ^(UpW, UaW) = UB a ^(p,cr)U^. Also, the value / does not change since the 
transformation affects both p and the basis that it is pinched onto in the very same manner. 

Normalization. Let d = 1. For every t £ [0,1] we have i? t (r||s) = — log(s/r), so that i? t (l||l/2) = 1 
and normalization is given. 

Order axiom. Let p < a. Then for every of the unitary matrices U t we get UtpUj < UtcrUj. 
Thus for every n £ N, t £ [0,1] as well as every pair (/, A), 

-^logi^Pf^p-^UtaUt)^} < -^logtviPf^UMp-^pUm (185) 

= -±logtr {PfAUtpUl) 9n }. (186) 

This implies that, for every t £ [0,1], 

Rt(p\\a)<D(p\\p) = 0 (187) 

and since D(p\\p) = 0 we get Rt(p\\cr) > 0. 

In case that p > a we get Rt(p\\cr) > 0 in the very same manner. 

Additivity. Is given by definition. 

Generalized mean value axiom. Since we consider the case d = 2 only, we get relieved from 
a heavy burden: Within B(C 2 ) the notion Rt(p © t||<t © w) implies that p , t, a and w are rank-one 
operators which satisfy p = r ■ |u)(u|, a = s|u)(|, r = t\v)(v\ and w = u>|v}(u| for thwo orthogonal and 
normalized vectors u,v £ C 2 . This immediately implies that [p + r, er + w] = 0 and also [p, a] = [r, w] = 0. 
This implies that all the functions Rt that occur are calculated as if they were classical Kullback-Leibler 
divergences. We thus see that the generalized mean value axiom cannot hold since it does not hold for 
the classical Kullback-Leibler divergence. 

Data Processing Inequality. Proving that DPI is valid (if that is true) seems a challenging 
and potentially fruitful task. While it is certainly clear that Ri = D satisfies the DPI, we are not yet 
aware of the methods which could be used to prove that DPI holds for other t £ [0,1]. One way to do so 
would certainly be to employ results from representation theory, while another obvious way would be to 
prove that R t equals D a z for certain choices of parameters. Be aware though that it has been proven in 
PP that DPI does not hold for t = 0. 

VII Conclusion 

We have brought forward our approach from |20| and proven that it leads to nontrivial connections 
between quantum information theory, representation theory and matrix analysis. Specifically, we have: 

1. delivered an operational interpretation for the limit lima-^i D a in Theorem Q] 

2 . and defined a new and nontrivial class {7?t}te[o,il of functions on qubits that are intimately con¬ 
nected to quantum relative entropies in Definition [9] Our ability to define these functions is based 
on Theorem [4] 

3. We used our approach to guess a nontrivial formula for minors of unitary matrices in Lemma [3j 

4. We have additionally provided explicit formulas for certain Hilbert-Schmidt scalar products in 
Theorems [2] and [3] These may turn out to be useful in later work. 
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We had to leave open the question of a more generic connection between the class {-Rt}te[o,i] and the 
a — z relative entropy as well as further connections to matrix geometry. It is our hope that such 
connections could lead to an expansion of our definition to arbitrary d and that this connection would in 
turn be able to provide meaningful statements on the intersection between quantum information theory 
and representation theory. 
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